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Spatial coherence resonance in a spatially extended system that is locally modeled by Hodgkin-Huxley (HH) 
neurons is studied in this paper. We focus on the ability of additive temporally and spatially uncorrelated 
Gaussian noise to extract a particular spatial frequency of excitatory waves in the medium, whereby examining 
also the impact of diffusive and small-world network topology determining the interactions amongst coupled 
HH neurons. We show that there exists an intermediate noise intensity that is able to extract a characteristic 
spatial frequency of the system in a resonant manner provided the latter is diffusively coupled, thus indicating 
the existence of spatial coherence resonance. However, as the diffusive topology of the medium is relaxed via 
the introduction of shortcut links introducing small-world properties amongst coupled HH neurons, the ability 
of additive Gaussian noise to evoke ordered excitatory waves deteriorates rather spectacularly, leading to the 
decoherence of the spatial dynamics and with it related absence of spatial coherence resonance. In particular, 
already a minute fraction of shortcut links suffices to substantially disrupt coherent pattern formation in the 
examined system. 

PACS numbers: 05.45.-a, 05.40.-a, 87.18.Sn, 89.75.Kd 
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Nontrivial effects of noise on nonlinear dynamics have 
been a vibrant topic for many years. It is thoroughly doc- 
umented and established that noise can play a construc- 
tive role in different types of nonlinear dynamical sys- 
tems. Stochastic and coherence resonance are just two, 
perhaps most prominent, examples of this fact. The no- 
tion of coherence resonance is particularly inspiring as it 
messages that an appropriate intensity of noise alone is 
sufficient to evoke ordered temporal responses of a non- 
linear dynamical system. Nowadays, however, effects of 
noise on spatially extended systems have gradually slipped 
into the focus of many scientists working in diverse fields 
of research, consequently spawning the need to investigate 
whether phenomena observed previously for isolated dy- 
namical system can also be observed, at least conceptually 
similar, if the latter are coupled. Indeed, it has been shown 
that the coherence resonance phenomenon originally re- 
ported for dynamical systems evolving only in time can 
also be observed in spatially extended systems that are 
locally described by excitable nonlinear dynamics. Im- 
portantly thereby is the fact that the previously studied 
order in the temporal dynamics has been "replaced" by 
the order of the noise-induced spatial dynamics, which 
mostly manifests as propagating waves of excitatory events 
throughout the spatial grid. Presently, we aim to extend 
the scope of spatial coherence resonance by confirming its 
possibility also in models of neuronal dynamics, in par- 
ticular by employing as the constitutive unit of the spa- 



tially extended system the renowned HH model. We show 
that while the diffusive connectivity of coupled neurons 
warrants the observation of noise-induced pattern forma- 
tion and with it related spatial coherence resonance, the 
small-world topology is not an appropriate medium for 
such observations. More precisely, even a minute fraction 
of shortcut links amongst distant neurons prohibits noise- 
induced waves to be ordered, and hence also precludes 
the observation of spatial coherence resonance. Since the 
present study is set-up around a comprehensive HH model 
of neuronal dynamics, the presented results should prove 
valuable not just from the purely theoretical point of view, 
but hopefully also from the experimental point of view, es- 
pecially by shedding light into the functioning of neural 
tissue. 



I. INTRODUCTION 

Randomness is a common feature in the real world. It is 
well known that noise can play a constructive role in differ- 
ent nonlinear dynamical systems, such as optical devices, 1 
electronic circuits, 2 or neural tissue. 3 Stochastic and coher- 
ence resonance are two prominent examples of such construc- 
tive effects that rely on the influence of noise on nonlinear 



systems. • Remarkably, by the coherence resonances^ an 
enhanced ordered behavior results solely from the introduc- 
tion of noise in the absence of additional weak determinis- 
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tic signals that are otherwise a standard ingredient by the ob- 
servation of stochastic resonance. 4 The phenomenon of co- 
herence resonance has been studied extensively also in non- 
identical coupled neurons 9 '- as well as one-dimensional ar- 
rays of nonlinear dynamical systems in genera h 11 ' 12 ' 13 Zhou 
and Kurths 14 showed, for example, that in a weakly coupled 
region noise could induce spatio-temporal coherence reso- 
nance but reduced the degree of synchronization by stronger 
couplings. 

Recently, many authors have shifted their interests to the in- 
fluences of noise on two-dimensional spatially extended sys- 
tems, and to their spatial rather than temporal dynamics in 
particular. 15 Spatio-temporal stochastic resonance has been 
reported first in, 16 while the spatial coherence resonance has 
also been studied in many different types of spatially extended 
systems, e.g. in chlorine dioxide-iodine-malonic reactions, 17 
Rulkov maps, 18 or excitable biochemical media. 19 More pre- 
cisely, Carrillo et al. 17 demonstrated that spatial coherence 
resonance could be evoked close to pattern-forming instabil- 
ities, mimicking one of the mechanisms of standard tempo- 
ral coherence resonance, while Perc et alJ&X? showed that in 
locally discrete excitable spatially extended systems and ex- 
citable biochemical media additive or internal noise could also 
lead to the observation of spatial coherence resonance. Im- 
portantly, however, there is still a lack of comprehensive stud- 
ies investigating whether the above results concerning spatial 
coherence resonance can be obtained also in more complex 
spatially extended systems, particularly such that are locally 
modeled by realistic nonlinear dynamical systems that faith- 
fully describe a real-life biological process, as is for example 
the neuronal activity. 

In this paper, we would like to address the above-described 
void and extend the scope of spatial coherence resonance 
by examining the possibilities of its existence in a spatially 
extended system that is locally modeled by excitable HH 
neurons. 20 More precisely, we analyze spatial frequency spec- 
tra of the examined medium in dependence on different lev- 
els of additive noise and topologies of the interaction network 
constituting the couplings amongst the HH neurons. By calcu- 
lating the average spatial structure function, we present con- 
clusive evidences for spatial coherence resonance provided 
the neurons are diffusively coupled. In particular, we show 
that then there exists an optimal level of additive noise for 
which a particular spatial frequency of excitatory events in 
the medium is best pronounced. We emphasize that thereby 
no additional deterministic inputs are introduced to the sys- 
tem, and the latter is locally initiated from steady state initial 
conditions. In contrast with results obtained when employ- 
ing diffusive nearest-neighbor interactions, however, we find 
that small-world interaction network s 21,22 fail to sustain co- 
herent patterns of spatial noise-induced excitations, and thus 
do not warrant the observation of spatial coherence resonance. 
More precisely, even a minute fraction of rewired links, con- 
necting two nearest neighbors in case of diffusive coupling, 
heavily disrupts coherent pattern formation in the medium and 
induces decoherence of excitatory waves. We thus reveal the 
impact of additive noise and small-world topology on an en- 
semble of HH neurons, hence providing interesting theoretical 



insights into to the functioning of neural tissue. 

The paper is structured as follows. Section II is devoted to 
the description of the employed spatially extended system and 
the HH mathematical model as its main ingredient. Sections 
III and IV feature evidences for spatial coherence resonance 
on diffusive grids and decoherence of pattern formation on 
small-world networks, respectively, while in the last Section 
we summarize the results and outline biological implications 
of our findings. 

II. MATHEMATICAL MODEL 

In 1952, Hodgkin and Huxley presented a mathematical 
model to predict the quantitative behavior of an isolated squid 
giant axon. 20 Since then, the model has become a paradigm 
for mathematically describing neuron functioning, and many 
authors have studied its nonlinear dynamics. In particular, it 
has been discovered that with the change of an external cur- 
rent injected into the cell membrane, different bifurcations 
can occur rendering a stable HH neuron oscillatory, bursting, 
or even chaotic. 23 ' 24 - 25 ' 26 ' 27 ! 28 Moreover, the dynamics of cou- 
pled HH neurons has also been studied quite extensively in 
the past i 29,30 ' 31 In this paper, we employ the HH model as the 
constitutive unit of a spatially extended system, and study the 
effects of different intensities of additive noise and topologies 
of the interaction network on the spatial dynamics of excita- 
tory waves. The equations of a single HH model are given as 
follows: 

CdV/dt = -g Na m 3 h(V - V Na ) - g K n\V - V K ) 

-9l(V - V L ) + I exU (la) 

dm/dt — a m (l — m) — [3 m m, (lb) 

dh/dt = a h (l-h)-0 h h, (lc) 

dn/dt — a n (l — n) — [3 n n. (Id) 

In Eq. ([TJ V is the membrane potential of the neuron, I ext is 
the external stimulus current, whereas m, h, and n are the gat- 
ing variables describing the ionic transport through the mem- 
brane. Moreover, the constants gpf a , gx, and g^ are maximal 
conductances for ion and leakage channels, while Vno,, Vk, 
and Vl are the corresponding reversal potentials. The func- 
tions a m , ah, ot n , (3 m , fth, and (3 n are defined as: 

a m = 0.1(V + 40)/{l - exp[-(V + 40)/10]}, (2a) 

a h = 0.07exp[-(V + 65)/20], (2b) 

a n =0.01(V + 55)/{l-exp[-(7 + 55)/10]}, (2c) 

P m = 4.0exp[-(V + 65)/18], (2d) 



P h = {l + exp[-(V + 35)/10}}-\ (2e) 

(3„ = 0.125 exp[-(V + 65)/80]. (2f) 

In the literature one may find different forms of the HH equa- 
tions, mainly depending on the zeros of the potential. Here 
we choose the form so that the membrane potential at rest 
is V « — 61mV. Then according to the studies of Hodgkin 
and Huxley, 20 the parameter values in Eq. (Q]i are as follows: 
C = 1/iF/cm 2 , g^a = 120ms/cm 2 , gx — 36ms/cm 2 , = 
0.3ms/cm 2 ,V Afa = 50mV, V K = -77mV, and V L = 
-54.4mV. 

First, we outline some properties of a single HH neuron 
evoked by the above parameters dependent on the external 
stimulus Iext^ When I ex t < Io ~ 6.2 there exists only a 
globally stable fixed point. The birth of stable and unstable 
limit cycles occurs at I ext — Io ~ 6.2 via a saddle-node bi- 
furcation. For Iq < I ext < 1% « 9.8 there coexist a stable 
fixed point, a stable limit cycle, and an unstable limit cycle. 
The unstable limit cycle constitutes the boundary separating 
the attractive basins corresponding to the fixed point and the 
limit cycle, respectively. For a more precise description of the 
bifurcation structure of the HH model, in dependence also on 
other system parameters, we refer to ) 23 ' 24,25 - 26 ' 27 ' 28 

Based on the original HH model given by Eqs. (fl]i and (O, 
we introduce the spatially extend system 

CdVij/dt = v„/»;. ( /', ,•: F,., - V Na ) 

-9Kr4,j( V i,j ~ V k) + f& ,j (t) 

-SX(K,j ~V L )+ Iext 

+D^e ilj M<y k ,i-V ilj ), (3a) 

k,l 

drrii.j/dt = a„ H ] (l - m it j) - P mitj mij, (3b) 
dh id /dt = a, H ] (l - h itj ) - Ph i:j hi,j, (3c) 

driij/dt = a ni j (1 - n it j) - f3 niij n id , (3d) 

where subscripts i, j = 1, . . . , N denote each of the N x N 
coupled HH neurons. The sum in Eq. (fjjl runs over all lat- 
tice sites whereby £i,j : k,i = 1 if the site (k, I) is coupled 
to and £ij,k,i = otherwise. If the fraction of ran- 

domly rewired links q, constituting the small-world topology 
if < q <C 1, equals zero, £i,j,k,l — 1 on ly if (k, I) is 
one of the four nearest neighbors of site on a regular 
two-dimensional mesh. Thereby, we obtain a diffusively cou- 
pled network of HH neurons as depicted in Fig. [TJ a )> which 
will form the backbone for the study of noise-induced spatial 
dynamics in Section III. If q > 0, however, the correspond- 
ing fraction of nearest-neighbor links is randomly rewired via 
the analogy with the creation of Watts-Strogatz small-world 
networks* 2 -!^ whereby we preserve the initial connectivity 
z = 4 (as by the diffusive nearest-neighbor coupling) of 
each HH neuron as depicted in Fig. [TJb) to focus explicitly 
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FIG. 1: Examples of considered network topologies. For clarity only 
a 6 x 6 excerpt of the whole network is presented in each panel, (a) 
Diffusively coupled network characterized by q — 0. Each vertex is 
directly connected only to its four nearest neighbors, hence having 
connectivity z = 4. (b) Realization of small-world topology via 
random rewiring of a certain fraction q of links, constrained only by 
the requirement that the initial connectivity z — 4 of each unit must 
be preserved. 

on the effect of network randomness rather than possible ef- 
fects originating from different numbers of inputs per neuron. 
The precise procedure for generating such small-world net- 
works is given in, 32 while results obtained when using them 
as interaction networks for the noise driven HH spatially ex- 
tended system will be presented in Section IV. Importantly, 
each interaction network was generated at the beginning of 
a particular simulation and was held fixed, and if necessary, 
results presented below were averaged over 30 different real- 
izations of the interaction network by each q. Turning back 
to Eq. (|3), cr is the main control parameter to be varied in 
this study, denoting the standard deviation of additive uncor- 
rected Gaussian noise £jj that satisfies < £i.j(t) >= and 

< ii, 3 {t)im,,n{t ) >= S(t-t')6i im Sj,n- Moreover, D = 0.35 
is the presently employed coupling strength, while all other 
parameters are taken the same as mentioned by the descrip- 
tion of a single HH neuron. In the present study the periodic 
boundary condition is used, namely: Voj = Vnj, Vn+ij — 
Vl,j> Vifi = Vi t N, and V^at+i = Vi t %, and the system size is 
N = 128 in each spatial dimension of the two-dimensional 
grid. In order to study explicitly the impact of noise on the 
spatial dynamics of the system, we take the external current 
Iext = 6.1 (in accordance with the above outlined properties 
of a single HH neuron assuring a unique stable steady state 
solution of each coupled unit) and initiate all units from the 
excitable steady state (-61.198, 0.08199, 0.46014, 0.37727). 
Thus, without the addition of noise the medium would remain 
quiescent forever, and consequently the dynamics of below- 
discussed spatial patterns is evoked solely by additive Gaus- 
sian noise. 



III. SPATIAL COHERENCE RESONANCE BY DIFFUSIVE 
COUPLING 

As announced, in this section we take q — whereby 
the coupling term DJ2 £ i.j,k,i(Vk,i — Vi.j) in Eq. © reads 




FIG. 2: Characteristic snapshots of the spatial profile of Vjj obtained for (a) a = 1.1, (b) a = 1.3 , and (c) a = 1.9. For the intermediate 
value of a the spatial dynamics is clearly most coherent, exhibiting ordered circular waves propagating through the spatial grid. 
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+ + Vij+i — 4Vij). We start lar average of the structure function according to the equation 



by visually examining three snapshots of the spatial profile 
of Vij obtained by three different values of a, as presented 
in Fig. |2] Evidently, there exists an intermediate value of a, 
at which nicely ordered circular excitatory waves propagate 
through the spatial domain [Fig. 12b)]. On the other hand, 
smaller or larger values of a clearly fail to have the same ef- 
fect evoking either only small-amplitude deviations from the 
steady state [Fig. a)] or rather violent and uncorrelated ex- 
citations throughout the spatial grid [Fig. Efc)], respectively. 
Hence, the visual inspection of the snapshots already gives 
some indication for a typical coherence resonance scenario 
for the spatial dynamics of the studied HH medium. 

To make the above observations quantitative, we calculate 
the spatial structure function of the Vi.j field according to 



P{k x ,ky) = (H 2 {k X ,ky)), 



(4) 



where H(k x ,k y ) is the spatial Fourier transform of the Vij 
field at a particular time t, and (...) is the ensemble average 
over noise realizations. The results from numerical simula- 
tions for different a are presented in Fig. [3] The three de- 
picted panels correspond to the same values of a as used al- 
ready in Fig. |2] Indeed, the results in Fig. [3] fully support 
our visual assessments, as it can be observed nicely that for 
small and large noise levels [Fig.|3ja) and (c), respectively] the 
presented spectra show no particularly expressed spatial fre- 
quency, although a close examination of P(k x , k y ) at a = 1.9 
still reveals some remnants of structure formation in the sys- 
tem. Only for intermediate levels of noise the spatial structure 
function develops several well-expressed circularly symmet- 
ric rings, indicating the existence of a preferred spatial fre- 
quency induced by additive Gaussian noise. As the noise level 
is increased, random fluctuations start to dominate the spatial 
dynamics and thus, similar as by small noise levels, the char- 
acteristic waterfall-like outlay of P(k x , k y ) vanishes and no 
preferred spatial frequency can be inferred. 

Next, we exploit the circular symmetry of P(k x , k y ) as pro- 
posed by Carrillo et al^ to obtain an estimate for the signal- 
to-noise ratio (SNR) of the spatial dynamics of excitatory 
events for different a. In particular, we calculate the circu- 



p(k) = / p(k)dn k , 



(5) 



where k — (k x ,k y ), and J\. is a circular shell of radius 
k = \k\. Figure a) shows p(k) for the three different a 
corresponding to the values used already in Figs. [2] and [3] 
The presented results reveal that there indeed exists a partic- 
ular spatial frequency, marked with the vertical dashed line 
at k = k max , that is resonantly enhanced for some interme- 
diate level of additive Gaussian noise. Note that subsequent 
local maxima are just higher harmonics of k max . To quan- 
tify the ability of each a to extract the characteristic spatial 
frequency of the medium more precisely, we calculate the 
SNR as the peak height at k max normalized with the back- 
ground fluctuations in the system; namely p(k max ) jp, where 
p = [p(k 

max 

Ak a ) + p{k 

max 

+ Afcf,)]/2 is an approxi- 
mation for the level of background fluctuations in the system, 
whereby Ak a and Afcb mark the estimated width of the peak 
around k max at the optimal a. This is the spatial counterpart 
of the measure frequently used for quantifying constructive ef- 
fects of noise in the temporal domain of dynamical systems, 4 
whereas a similar measure for quantifying effects of noise on 
the spatial dynamics of spatially extended systems was also 
used in. 17 FigureUb) shows how the SNR varies with a. It is 
evident that there exists an optimal level of additive Gaussian 
noise for which the peak of the circularly averaged structure 
function is best resolved, thus clearly indicating the existence 
of spatial coherence resonance in the diffusively coupled HH 
medium. 

We argue that the existence of spatial coherence resonance 
in the studied HH medium, and with it related remarkable or- 
der of the spatial dynamics of noise-induced excitatory events 
observed by intermediate a, must be primarily attributed to 
the characteristic noise-robust excursion time inherent to all 
excitable systems. 8 The noise-robust excursion time, together 
with the rate of diffusive spread that is proportional to vD, 
introduces an eigenfrequency of spatial waves that can be res- 
onantly enhanced by an appropriate intensity of noise, as re- 
cently shown analytically for a simple toy model of excitable 
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FIG. 4: Evidences for spatial coherence resonance in the studied HH medium, (a) The circular averages of structure functions presented in 
Fig. [3] The dashed vertical line at k — k max marks the spatial frequency of excitatory events that is resonantly enhanced by an intermediate 
level of additive Gaussian noise, (b) The signal-to-noise (SNR) ratio in dependence on a (the curve is just a guide to the eye). 



dynamics. Our results thus extend the scope of spatial co- 
herence resonance by showing that real-life motivated more 
complex models of excitable dynamics, such as the HH math- 
ematical model, 20 fully conform to findings obtained on sim- 
pler models, and that indeed this phenomenon may have im- 
portant implications in several biological systems. 



IV. DECOHERENCE OF SPATIAL DYNAMICS BY 
SMALL- WORLD INTERACTIONS 

Next, it is of interest to investigate what are the impacts 
of small-world topology on the above-described phenomenon 
of spatial coherence resonance in the diffusively coupled HH 
medium. For this purpose we set q > and examine three 
characteristic snapshots of Vij obtained by a = 1.3, which 
according to the results presented in Section III is the optimal 
noise intensity warranting superbly ordered spatial dynamics 
of excitatory events. Noteworthy, since we consider regular 



small-world networks warranting an equal number of inputs 
(z = 4) to each coupled neuron, the optimal a for coherent 
pattern formation (if at all possible) does not depend on q, as 
will be confirmed also in Fig. |6ja) below. The snapshots pre- 
sented in Fig. [5] clearly evidence that increasing values of q 
progressively hinder coherent pattern formation. In particular, 
while for q = 0.0001 [Fig.|5ja)] the spatial periodicity of ex- 
citatory waves is still eligible, the same is much less true for 
q = 0.0008 [Fig.|5jb)], and is certainly completely false for 
q = 0.002 [Fig. He)]. Thus, if as little as 0.08 % of nearest- 
neighbor links are rewired the pattern formation is impaired, 
while by w 0.2 % noisy perturbations are completely unable 
to induce coherent spatial dynamics of excitatory events. In- 
deed, the snapshot presented in Fig. |3c) lacks any ordered 
structure, although the same intensity of noise applied to a 
diffusively coupled HH medium evokes superbly ordered cir- 
cular waves [see Fig. |2jb)] evidencing fascinating examples 
of spatial pattern formation out of noise. 

To give a quantitative analysis of the observed phenomenon 



(a) 



(b) 






FIG. 5: Characteristic snapshots of the spatial profile of Vi ,j obtained at the optimal value of a — 1.3 for (a) q — 0.0001, (b) q — 0.0008, and 
(c) q = 0.002. 



outlined in Fig. [5] we calculate the circularly average struc- 
ture function p(k) of the spatial dynamics in dependence on 
a and q. The inset in Fig. |6jb) shows three p(k) obtained by 
the same parameter values yielding snapshots of spatial pro- 
files shown in Fig. [5] In accordance with the interpretation of 
results presented in Section III, it can be inferred that increas- 
ing values of q indeed impair the ability of noise to evoke or- 
dered patterns with a predominant spatial frequency, because 
the peak height of p(k) at k = k. max becomes smaller in com- 
parison to the level of background noise as the fraction of 
rewired links increases. While by q = 0.0001 and q = 0.0008 
the predominant spatial periodicity is still present, it vanishes 
completely by q = 0.002 as the outlay of p(k) becomes essen- 
tially flat. These observations can be captured succinctly by 
calculating the SNR in dependence on a by different q. The 
results presented in Fig. |6ja) evidence clearly that the phe- 
nomenon of spatial coherence resonance fades continuously 
as q increases, disappearing almost completely by q = 0.002. 
As already noted, however, the optimal a for coherent pattern 
formation (as much as allowed by the topology of the under- 
lying interaction network) remains virtually identical by all q, 
equaling a = 1.3, as denoted by the dashed vertical line in 
Fig. |6ja). We exploit this fact in Fig. |6jb), where the SNR 
is plotted in dependence on q for the optimal a. Clearly, the 
SNR decreases continuously as q increases, thus evidencing 
the decoherence of noise-induced spatial dynamics due to the 
introduction of small-world topology in the studied HH spa- 
tially extended system. 

While the existence of spatial coherence resonance in the 
diffusively coupled HH medium was attributed to the char- 
acteristic noise-robust excursion time of the local excitable 
dynamics, 8 which together with the rate of diffusive spreading 
proportional to \[D introduces an eigenfrequency, or charac- 
teristic spatial scale, of waves that can be resonantly enhanced 
by an appropriate intensity of noise, 33 the decoherence of the 
spatial dynamics due to the introduction of small-world topol- 
ogy is simply a consequence of the disruption of this inherent 
spatial scale. 34 In particular, while the excursion time of in- 
dividual space units remains unaltered by the introduction of 
long-range couplings, the spread rate of excitations is indi- 
rectly very much affected by increasing values of q. Namely, 



the introduction of shortcut links decreases the typical path 
length between two arbitrary sites in comparison to a diffu- 
sively coupled network. Thus, while the spread rate of exci- 
tations is still proportional to y/D, the typical path length be- 
tween two arbitrary grid units decreases dramatically, which 
in turn has the same effect as if D would increase. Therefore, 
in a small-world network, a locally induced excitation can in- 
stantly reach much more distant sites than in case of strict 
diffusive coupling, in turn facilitating noise-induced synchro- 
nization of distinct network units,— or assuring fast response 
abilities of the system.- 36 Importantly, however, the typical 
path length between two arbitrary sites decreases only on av- 
erage, meaning that there does not exist an exact path length 
defining the distance between all possible pairs of sites in a 
small-world network. Thus, a given local excitation can, dur- 
ing the excursion time, propagate to the most distant site or 
just to its nearest neighbor, whereby clearly the well-defined 
inherent spatial scale existing for the diffusive coupling is 
lacking, and this leads to spatial decoherence of noise-induced 
spatial patterns already for very small values of q, as empha- 
sized throughout this Section. 



V. SUMMARY 

We have shown that temporally and spatially uncorrelated 
additive Gaussian noise is able to extract a characteristic spa- 
tial frequency of an excitable HH medium in a resonant man- 
ner if the latter is coupled diffusively. In particular, it then 
exists an optimal level of additive noise for which the spatial 
dynamics of the system is highly coherent and superbly or- 
dered. Thereby, no additional deterministic inputs were intro- 
duced to the system and all units were initiated from steady- 
state excitable conditions. Thus, the results presented in Sec- 
tion III offer convincing evidences for the existence of spatial 
coherence resonance in the diffusively coupled HH medium. 
On the other hand, when the interactions amongst coupled 
neurons are governed by a small-world network, the ability 
of noise to induce ordered spatial dynamics vanishes quickly 
as q increases. Our results suggest that as little as 0.2 % 
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FIG. 6: Evidences for decoherence of the spatial dynamics in the studied HH medium with small- world topology, (a) The signal-to-noise 
ratio (SNR) in dependence on a for different q. (b) SNR in dependence on q by the optimal a denoted by a dashed vertical line in panel (a). 
Inset features the circular average of the structure function p(k) for the optimal a and three different values of q, corresponding to the three 
snapshots presented in Fig. [5] The values depicted in the main panels resulted from averaging the results over 30 different realizations of the 
interaction network by each q, whereas curves are just guides to the eye. 



of rewired nearest-neighbor links suffice to completely hin- 
der coherent pattern formation out of noise within the cur- 
rently employed setup. Noteworthy, we have performed addi- 
tional simulations using the FitzHugh-Nagumo system^! and 
the Goldbeter model for calcium oscillations 38 as governing 
units of the medium, whereby the former excitable dynam- 
ics is governed by a steady node while the latter is charac- 
terized by a stable focus situated prior to a Hopf bifurcation. 
Irrespective of these particularities, we were able to observe 
coherent patterns emerging out of noise provided the units 
were diffusively coupled. Conversely, if small-world connec- 
tivity was introduced the ordered waves vanished quickly al- 
ready by q <C 1, albeit stable foci could sustain some order 
in the spatial dynamics up to 0.5 % of rewired links due to 
larger re-settlement times associated with their excitable dy- 
namics. Larger re-settlement times following an excitation 
yield waves with larger wavelengths, and hence then the prob- 
ability of shortcuts to link two quiescent units increases, thus 
leaving the spatial dynamics temporarily unaffected. Aside 
from these rather minute differences, however, we conclude 
that the present findings appear to be robust with respect to 
variations in models governing the local excitable dynamics 
of diffusive and small-world networks. 

Since the present results rely on the use of the HH model for 
the description of the local dynamics of the excitable medium, 



the present study might also have some biological implica- 
tions. Specifically for neural systems, it has been argued that 
excitable media guarantee robust signal propagation through 
the tissue in a substantially noisy environment. 39 It would thus 
be fascinating to elucidate if spatial coherence resonance in 
the neural system can be confirmed also experimentally. The 
above theoretical results indicate that such findings might in- 
deed be attainable at least if the units are diffusively coupled, 
and that it thus seems reasonable to pursue this problem in 
the future. Moreover, our results obtained when using small- 
world networks may help by revealing structural functionality 
of complex topologies within the neural apparatus, although 
additional studies regarding the concept of function-follow- 
form 4 ' 41 ! 42 ' 43 in the presence of noise and other uncertainties 
are necessary to clarify the importance of different structural 
physiological properties of such networks. 



Acknowledgments 

X. Sun and Q. Lu are grateful for the support of the Na- 
tional Science Foundation of China (Grants 10432010 and 
10702023). M. Perc acknowledges support from the Slove- 
nian Research Agency (Grant Zl-9629). 



* Electronic address: sxiaojuan@ss.buaa.edu.cn 

' Electronic address (corr.): matjaz.perc@uni-mb.si 

* Electronic address: qishaolu@hotmail.com 

§ Electronic address: jkurths@agnld.uni-potsdam.de 

1 A. Fioretti, L. Guidoni, R. Mannella, and E. Arimondo, J. Stat. 
Phys. 70, 403 (1993). 

2 T. Zhou, F. Moss, and P. Jung, Phys. Rev. A 42, 3161 (1990). 



3 E. Ullner, A. Zaikin, J. Garcfa-Ojalvo, and J. Kurths, Phys. Rev. 
Lett. 91, 180601 (2003). 

4 L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni, Rev. Mod. 
Phys. 70, 223 (1998). 

5 B. Lindner, J. Garcfa-Ojalvo, A. Neiman, and L. Schimansky- 
Geier, Phys. Rep. 392, 321 (2004). 

6 G. Hu, T. Ditzinger, C.Z. Ning, and H. Haken, Phys. Rev. Lett. 



8 



71, 807 (1993). 

7 W.J. Rappel and S.H. Strogatz, Phys. Rev. E 50, 3249 (1994). 

8 A.S. Pikovsky and J. Kurths, Phys. Rev. Lett. 78, 775 (1997). 

9 C.S. Zhou, J. Kurths, and B. Hu, Phys. Rev. Lett. 87, 098101 
(2001). 

C.S. Zhou and J. Kurths, Chaos 13, 401 (2003). 

1 H.S. Wio, Phys. Rev. E 54, R3075 (1995). 

2 S.K. Han, T.G. Yim, D.E. Postnov, and O.V. Sosnovtseva, Phys. 
Rev. Lett. 83, 1771 (1999). 

3 A. Neiman, L. Schimansky-Geier, A. Cornell-Bell, and F. Moss, 
Phys. Rev. Lett. 83, 4896 (1999). 

4 C.S. Zhou and J. Kurths, Phys. Rev. E 65, 040101 (R) (2002). 

5 F. Sagues, J.M. Sancho, and J. Garcfa-Ojalvo, Rev. Mod. Phys. 
79, 829 (2007). 

6 P. Jung and G. Mayer-Kress, Phys. Rev. Lett. 74, 2130 (1995). 

7 O. Carrillo, MA. Santos, J. Garcia-Ojalvo, and J.M. Sancho, Eu- 
rophys. Lett. 65, 452 (2004). 

8 M. Perc, Chaos, Solitons and Fractals 31, 64 (2007). 
M. Gosak, M. Marhl, and M. Perc, Biophys. Chem. 128, 210 
(2007). 

20 A.L. Hodgkin and A.F Huxley, J. Physiol. 117, 500 (1952). 

21 D.J. Watts and S.H. Strogatz, Nature 393, 440 (1998). 

22 D.J. Watts, Small world: The Dyanmics of Networks Between 
Order and Randomness (Princeton University Press, Princeton, 
1999). 

23 S.G Lee, A. Neiman, and S. Kim, Phys. Rev. E 57, 3292 (1998). 

24 J. Rinzel, Fed. Proc. 37, 2793 (1978). 

25 W.C. Troy, Q. Appl. Math. 36, 73 (1978). 

26 B. Hassard, J. Theory Biol. 71, 401 (1978). 



27 H. Fukai, T. Nomura, S. Doi, and S. Sato, Biol. Cybern. 82, 223 
(2000). 

28 J. Wang, J.M. Geng, and X.Y. Fei, Chaos, Solitons and Fractals 
23, 973 (2005). 

29 D. Hansel, G. Mato, and C. Meunier, Europhys. Lett. 23, 367 
(1993). 

30 M.S. Wang, Z.H. Hou, and H.W. Xin, Phys. Lett. A 334, 93 

(2005) . 

31 Q.Y. Wang, Q.S. Lu, and G.R. Chen, Eur. Phys. J. B 54, 255 

(2006) . 

32 G. Szabo, A. Szolnoki, and R. Izsak, J. Phys. A: Math. Gen. 37, 
2599 (2004). 

33 M. Perc and M. Marhl, Phys. Rev. E 73, 066205 (2006). 

34 M. Perc, New J. Phys. 7, 252 (2005). 

35 Z. Gao, B. Hu, and G. Hu, Phys. Rev. E 65, 016209 (2001). 

36 L.F. Lago-Fernandez, R. Huerta, F. Corbacho, and J.A. Sigiienza, 
Phys. Rev. Lett. 84, 2758 (2000). 

37 R. Fitzhugh, Biophys. J. 1, 445 (1961). 

38 A. Goldbeter, G. Dupont, and M.J. Berridge, Proc. Natl. Acad. 
Sci. USA 87, 1461 (1990). 

39 E.M. Izhikevich, Int. J. Bifurcation Chaos 10, 1171 (2000). 

40 V.L. Towle, R.K. Carder, L. Khorasanil, and D. Lindberg, J. Clin. 
Neurophysiol. 16, 528 (1999). 

41 C.C. Hilgetag, GA. Burns, MA. O'Neill, J.W. Scanell, and 
M.P. Young, Phil. Trans. R. Soc. B. 355, 91 (2000). 

42 R. Segev, M. Benveniste, Y. Shapira, and E. Ben- Jacob, Phys. Rev. 
Lett. 90, 168101 (2003). 

43 V. Volman, I. Baruchi, and E. Ben- Jacob, Phys. Biol. 2, 98 (2005). 



